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An instructional laboratory experiment to measure the muon mass is described. Using 
coincidence-anticoincidence detection, the decay of a cosmic-ray muon into an electron (or positron) 
is observed in a multiplate spark chamber, and recorded with a triggered CCD detector. The energy 
of the charged decay-product particle is then quantified by counting the number of spark gaps ns it 
traverses before being stopped by the chamber's aluminum plates. By running this apparatus under 
computer-control for several hours, the number of product-particles N(ns) with various n s -values 
is obtained. The muon mass is obtained by a least-squares fit, which compares the experimentally 
observed N(ns) with simulation values predicted for this distribution by the Fermi description of 
muon decay via the weak interaction. We present the results for the muon mass we have obtained 
and discuss the simulation we developed to account for the observed skewing of N(ns) due to the 
various directions the spark-producing product particles move as well as the escape of some of the 
higher-energy particles from the chamber. 



I. INTRODUCTION 

Earth is continually bombarded by cosmic rays, high- 
speed subatomic particles produced by astrophysical pro- 
cesses such as supernovae. Many of the first experiments 
in elementary particle physics used cosmic rays as their 
particle source, and even today, they continue to be stud- 
ied in experiments that probe the very highest particle 
energies. While this source of high-speed particles is om- 
nipresent and available at no cost, the flux of cosmic rays 
is relatively low. Thus experiments done with cosmic rays 
can best probe processes with high intrinsic rates. 

For over forty years, the presence of high-speed muons 
at Earth's surface has allowed educators to design in- 
structional laboratory experiments that illustrate the 
concerns and experimental methods of elementary parti- 
cle physics. These muons are the byproduct of collisions 
between cosmic rays and nuclei near the top of Earth's 
atmosphere. By far, the most widely adopted such ex- 
periment measures the muon's lifetime [TJ [2j 02 HJ [5j |6] , 
but others include observation of the time dilation ef- 
fect |6j [7] and determination of the muon's magnetic 
moment. [8 In 1964, researchers at the University of 
Michigan suggested that the mass of the muon could be 
measured through an experiment that analyzed tracks 
produced by muon decay within a spark chamber. [9] 
These researchers built a prototype system and demon- 
strated the feasibility of this approach. However, at that 
time, the necessary setup required too much specialized 
technical expertise for most educators to copy and the ex- 
periment was too cumbersome to perform because track 
images were captured on photographic film. Thus, this 
experiment has never been embraced by instructional lab 
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developers. In this paper, we describe the muon mass 
experiment and how it can be implemented today with 
relative ease through the use of modern instrumentation, 
including CCD-based image acquisition. We go on to 
show how an accurate value for the muon mass can be 
determined from the data gleaned from this experiment 
using an easy to replicate software simulation program. 

II. THEORY 

A muon is a negatively charged elementary particle 
from the lepton family with a mass that is about 200 
times that of an electron. With effectively 100% prob- 
ability, a free muon decays into an electron, an electron 
antineutrino, and a muon neutrino with a mean lifetime 
of 2.2 /is. That is, the muon's principal decay mode is 
ji~ — > e~ + V e + v^. Similarly, the principal decay mode 
of the muon's antiparticle is /i + — > e + + v e + V^. As 
a result of the interaction of cosmic rays with the nu- 
clei of air molecules at the top of Earth's atmosphere, 
Earth's surface is bombarded by a stream of high-speed 
charged secondary cosmic ray particles, 75% of which are 
a roughly equal mixture of muons and antimuons. In the 
vertical direction, the total flux per unit solid angle of 
these secondary particles is about 0.66/cm 2 • sr • min. pQ 
Using this value, along with the empirically determined 
cosine-squared variation of particle flux with angle from 
the vertical, it can be shown that the flux of charged sec- 
ondary cosmic-ray particles impinging on a horizontally 
aligned detector is on the order of 1/cm 2 • min. 

The goal of our experiment is to stop a significant 
number of secondary cosmic-ray muons within a spark 
chamber and observe each decay into an electron (to save 
words, we will simply describe the decay of a muon into 
an electron, but our discussion likewise applies to the 
decay of an antimuon decay into a positron). The spark 
chamber is a charged-particle detector that makes visible 
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the muon's path prior to being stopped as well as the path 
of the product electron after the decay. Two product 
neutrinos are also produced by the muon's decay, but be- 
cause they are neutral their paths are not recorded by the 
spark chamber. Our chamber consists of a vertical stack 
of aluminum plates. Most cosmic-ray muons that enter 
the chamber have energies greater than 1 GeV and so, al- 
though slowed by the plates, will pass completely through 
the chamber. However, a small population of low-energy 
muons will be effectively "stopped" (i.e., slowed to ki- 
netic energies on the order of the aluminum's ionization 
energy) as they pass through the sequence of plates, sup- 
pressing the time-dilation effect that allowed these short- 
lived particles to traverse the height of the atmosphere. 
Roughly half of these "stopped" muons then decay freely 
via fi~ — > e~ +V e + while the other half are captured 
by atoms and undergo the process of inverse beta decay 
fi~ + Al — > e~ + + Mg* (which is undetected by the 
spark chamber). On the other hand, all of the stopped 
antimuons decay freely via /i + — > e + + v e + V^. 

Consider the decay of a stationary muon of mass m^c 2 
(consistent with the shorthand used in the field of parti- 
cle physics, we will call this quantity the "muon mass," 
although it is actually the "muon rest energy") into an 
electron, an electron antineutrino, and a muon neutrino. 
Applying conservation of energy to this process, we find 
m M c 2 = E e + Ey e + E v , where the energies E ei E„ e1 
and E v of the three product particles can be taken (to 
a good approximation) to be entirely kinetic energy be- 
cause each of these particles has a rest mass much less 
than m^c 2 . Additionally, applying conservation of mo- 
mentum, we get = p e +p„ e +Vv^, where p e , p„ e , and 
p u are the relativistic momenta of the product parti- 
cles. With three product particles, the kinetic energy of 
the electron is not uniquely determined by the conser- 
vation laws, but instead a range of £^ e -values is possi- 
ble. The minimum value of this range is E e = (ne- 
glecting the electron's mass) in which case the electron 
is at rest and the two neutrinos move oppositely directed 
with equal momentum magnitude. At the other extreme, 
the maximum electron energy occurs when the two neu- 
trinos both move in the opposite direction to the elec- 
tron. Then, assuming the product particles are all mass- 
less, it is easy to show from the conservation laws that 
E e = \m^c 2 Thus, the allowed range of electron ener- 
gies is < E e < ^m^c 2 . The Fermi description of muon 
decay via the weak interaction, [11] in which the decay 
(in lowest order) is mediated by a W boson of mass M w , 
yields the following relation for the probability per unit 
time dT that a product electron will be produced with 
energy in the interval E e to E e + dE e 

dE e = {j^) 384^ (^^ E ^j ( 3 " 4 ^) 

< E e < X -m^c 2 (1) 
where g w is the weak coupling constant. Higher order 
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FIG. 1: Experimental set-up, which includes multiplate spark 
chamber, coincidence/anticoincidence electronics to detect 
muon decay within chamber, fast high- volt age pulser to cre- 
ate spark trail, and triggered CCD-based image acquisition 
to store data. 



electromagnetic effects modify this expression slightly 
(by approximately 1%). Noting that dT/dE e is inde- 
pendent of angle, Eq. predicts that, within a given 
time interval and solid angle, the probability P(E e ) of 
product electrons being produced with energy E e is 

P(E e ) = C(m M c 2 £ e ) 2 (3 - 4£ e /m M c 2 ) (2) 

where C is a constant. Hence, by fitting an experimen- 
tally determined plot of P(E e ) vs. E e to this predicted 
functional form, a value for m^c 2 can be determined. 



III. SYSTEM OVERVIEW 

Our system consists of four main components: spark 
chamber, particle-decay detection electronics, high- 
voltage pulsing electronics, and CCD-based image acqui- 
sition (see Fig. [I]). A short description of each of these 
components follows; more detailed descriptions are given 
in the Appendix. 

The spark chamber consists of a stack of 21 disk- 
shaped aluminum plates. Each plate has thickness t p = 
0.9525 cm (= 3/8 inch) and radius R = 7.62 cm (= 
3 inch) and is separated from its neighboring plates by 
"spark gaps" of thickness t g = 0.635 cm (= 1/4 inch) 
filled with a noble gas (in our case, neon). The odd- 
numbered plates are electrically grounded, while the 
even-numbered plates are connected to a fast high- 
voltage, high-current pulser. When a charged particle 
moves within the chamber, it ionizes the noble gas along 
its path, leaving behind a trail of ionized atoms and lib- 
erated electrons. If the pulser can be triggered to apply 
sufficiently high voltage (typically 6000 V) to the even- 
numbered plates more quickly than the mean ion-electron 
recombination time (on order of 5 /is), sparks will form 
in the gaps between the plates along the low-resistance 
ionized track, making visible the trajectory of the moving 
charged particle. 
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The decay of a cosmic-ray muon within the spark 
chamber is detected through the use of three scintillation 
detectors stacked vertically about the chamber. Each 
detector consists of a horizontally aligned light-tight 20 
cm x 20 cm x 2 cm plastic scintillator coupled to a 
blue-light sensitive photomultiplier tube (PMT). When a 
high-speed charged muon passes through such a detector, 
it causes the scintillator material to fluoresce [10] and the 
resulting photons are detected and converted to a short 
(« 5 ns), negative-going electrical signal by the attached 
PMT. This signal is amplified and then passed through a 
discriminator which converts the PMT signal into a short 
digital pulse. Two of these detectors (henceforth called 
X and Y) are placed above the chamber, separated by 
a vertical distance of approximately 1 meter, while the 
third detector (Z) is placed slightly below the chamber. 
When a muon with (close to) normal incidence enters 
the chamber from above, and then subsequently decays 
within the chamber, digital pulses will be produced by X 
and Y, while a digital pulse will not be produced by Z, 
assuming the decay-product electron is stopped within 
the chamber. Thus, if the digital outputs of the three 
discriminators are input to a logic unit that performs the 
coincidence/anticoincidence operation X • Y • Z (X AND 
Y AND notZ), the output of this logic unit signals the 
union's decay within the spark chamber. 

The X Y Z output of the logic unit is connected 
to the input of a high-speed pulser ("driver"). This 
driver unit produces the proper voltage pulse to trigger 
a high- voltage, high-current switch ( "thyratron" ) , which 
then causes a bank of capacitors to discharge, applying 
a (negative) pulse of several thousand Volts to the even- 
numbered plates of the spark chamber. 

Finally, a video camera focused on the spark chamber 
is configured to send its analog video signal to the video 
input of a frame grabber board within the expansion slot 
of a PC, while the X Y Z output of the logic unit is con- 
nected to the board's trigger input. The frame grabber is 
controlled via a Lab VIEW program, which acquires and 
stores a single image of the chamber whenever the logic 
unit detects a muon decay event. 



IV. EXPERIMENTAL PROCEDURE 

With high voltage applied to all three scintillation de- 
tector PMTs, each PMT output (after passing through 
a preamplifier) is input to a separate channel of a quad 
constant fraction discriminator (CFD). Using an oscil- 
loscope, the width of NIM digital pulses output by 
each CFD channel (in response to PMT input signal 
pulses) is adjusted to be 300 ns. In our setup, each 
detector's cross-sectional area A = 20 cm x 20 cm = 
400 cm 2 , so the rate of digital pulses produced by sec- 
ondary cosmic ray particles on each channel is about 
(1/cm 2 • min)(400 cm 2 )(lmin/60 s) w 7/s. 

The digital-pulse outputs of the CFD's three chan- 
nels X, y, and Z are input into the logic unit. The 



TTL output of this logic unit is connected to the high- 
speed driver, which triggers the spark chamber opera- 
tion. Proper functioning of the setup can be verified 
by programming the logic unit to perform the logic op- 
eration X • y, so that the spark chamber is triggered 
for any charged particle (i.e., not just those that decay) 
that passes into the chamber through the top two de- 
tectors. The expected rate of such X ■ Y coincidence is 
found as follows: X and Y each have area A = 400 cm 2 , 
and are separated by distance d = 100 cm. Thus, at 
each point on Y, the solid angle subtended by X is 
~ (loo cm) 2 = 0-04 sr, and so the expected rate at 
which vertically downward secondary cosmic ray parti- 
cles (75 % of which are muons and antimuons) will pass 
through both detectors producing coincidence signals is 
(0.66/cm 2 • sr • min)(400 cm 2 ) (0.04 sr) = 10/min. 

To execute a data run, the logic unit is configured to 
perform the logic operation X • Y • Z, the video camera 
is focused on the spark chamber along with its reflection 
in a nearby mirror angled at 45° (to obtain an orthog- 
onal viewpoint of the spark tracks), and a LabVIEW- 
based image acquisition program is started. Each time 
the logic unit triggers the spark chamber, it also sends a 
TTL trigger pulse that triggers the frame grabber board 
to acquire the video camera's current frame and save it 
as a jpeg file on the computer's hard drive with a unique 
file name. The experiment is run under computer con- 
trol in a darkened room overnight, typically resulting in 
about 100 saved images per hour. 

After a data run, all the images are inspected. Fig. 
[2] shows a sample image of a muon decay event, which 
consists of a direct and an orthogonal view of the cham- 
ber. The orthogonal view is the image produced by a 
45°-angled mirror located at the left side of the chamber; 
this image is smaller than the direct view because it is 
farther from the camera than the chamber. In Fig. [2j 
the muon enters from the top of the chamber and, after 
it decays, the product-electron moves downward to the 
left and toward the back of the chamber. To analyze this 
decay image, we construct two lines, one delineating the 
path of the incoming muon, and the other marking the 
path of the product electron. The kink at the intersec- 
tion of these lines defines the location at which the muon 
decayed. Once the point of decay is established, one then 
counts the number of sparks ns created by the electron 
as it is brought to rest (or else escapes from the chamber) 
during its traverse through the aluminum plates. In Fig. 
n s = 6. 

In contrast to the image in Fig. |2j many of our ac- 
quired images do not clearly show a muon decay event 
for the following reasons. First, although from Eq. ([!]), 
we expect an isotropic distribution of product-electron 
tracks, i.e., as many images with an electron going up 
as down from the decay site, we never get a clear decay 
image showing an upward moving electron. This observa- 
tion is simply an artifact of our spark-chamber detection 
method; to record a decay event with an upward moving 
electron would require each spark gap above the decay 
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• Experiment 
□ Simulation 



m^c = 100 MeV 
9 =30° 



012345678 
Number of Sparks 

FIG. 3: Comparison of experimental data (Nl°p — 43) with 
the best-fit result of the simulation, assuming m^c 2 = 100 
MeV and Omax — 30°. Error bars on the data are computed 
as the square root of the associated count number. 




the number of muon decays N exp (ns) in which product 
electrons produce ns sparks, and then plot N exp (ns) vs. 
ns as in Fig. [3] 



FIG. 2: Typical spark image of a muon decay with the direct 
and orthogonal (created by 45°-angled mirror) view of the 
chamber on right and left, respectively. Here, the muon enters 
from the top and, after it decays, the product electron moves 
leftward and toward the back of the chamber while producing 
6 sparks. Superimposed on each view is a software-generated 
outer and inner rectangle defining the boundaries of the entire 
chamber and the half-radius volume, respectively. The upper 
portion of the half-radius volume is the fiducial volume within 
which a valid decay must initiate. 



site to produce two sparks, one for the downward inci- 
dent muon and another for the upward directed electron. 
Unfortunately, each spark gap can only produce a single 
spark. Second, for a downward directed electron, we find 
that a spark is created only if the electron's path devi- 
ates up to roughly 30° from the chamber's axis. This 
observation results from the fact that the electric field 
produced by the plates is along the chamber's axis (i.e., 
normal to the plates) and if the ionized trail is inclined 
too steeply from this direction, no spark forms. Thus, in 
our inspection process, we select out the images which 
clearly display muon decay events. Because of the de- 
tails of spark formation described above, these images 
always have downward directed product-electrons, mov- 
ing at some angle less than about 30° to the chamber's 
axis. 

For reasons delineated in the next section, we analyze 
only the subset of those images for which the muon decay 
occurs within the chamber's "fiducial volume," which we 
define to be a cylinder in the top half of the chamber 
with a radius Rfid = R/2 measured from the chamber's 
central axis. After this selection process, one is left with 
a total of Nl° P images that record muon decays (typi- 
cally iVg°* « 50). Using these data, we then determine 



V. ANALYSIS OF DATA 

The product electron's initial kinetic energy E e (i.e., 
its energy immediately after the muon decay) can be de- 
termined by analyzing its path through the chamber's 
sequence of aluminum plates and spark gaps. Assum- 
ing a ballpark value of m M c 2 « 100 MeV, Eq. 
indicates that E e will be in the range from to 50 
MeV. As electrons with initial energies in this range are 
slowed during their travel through aluminum, the two 
most significant energy-loss processes are ionization and 
bremsst ranking, [12] whose stopping powers S = —dE/dl, 
i.e., energy loss per unit length, are approximately a con- 
stant and proportional to instantaneous energy E, re- 
spectively. In particular, the total stopping power for an 
electron in aluminum [I3j [14] is given by 
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where S = 5.09 MeV/cm and X Q = 8.9 cm is the "radi- 
ation length" for aluminum. Note that Si on iz = Sbrem at 
the "critical energy" E crit = S X = 45 MeV, and that 
for electron energies much less (much greater) than this 
critical energy, ionization (bremsstrahlung) is the domi- 
nant energy- loss process. Since the electron energies in 
our experiment do not deviate greatly from E critl we 
must account for both energy- loss mechanisms; thank- 
fully, it is easy to show from Eq. Q that the path length 
I in aluminum required to bring an electron with initial 
energy E e to rest is 
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The ideal experimental procedure would be as follows: 
For each observed muon decay, measure the length I of 
the product electron's path through the chamber's alu- 
minum plates, and then use Eq. Q to calculate the elec- 
tron's initial energy E e . Repeating this process for all 
the observed decays, determine the probability P(E e ) of 
a product particle being produced with energy E e , then 
fit the resulting plot of P(E e ) vs. E e to Eq. ([2]), and ex- 
tract a value for the muon mass from the fit. Because we 
obtain two perpendicular views of the electron's path, in 
principle, we could carry out this ideal procedure. How- 
ever, in practice, it proves very difficult to determine the 
length of the electron's path in three dimensions from our 
data, and so we take another tack. 

As mentioned in the previous section, we simply char- 
acterize the path of a product electron by counting the 
number of sparks ns it creates while moving in the cham- 
ber. Clearly, there is not a one-to-one correspondence be- 
tween ns and E e , e.g., with the same initial energy E e , 
an electron passing through the aluminum plates at nor- 
mal incidence will produce more sparks than one passing 
through the plates at a large incident angle (because of 
the increased path length per plate). However, by writ- 
ing a straightforward Monte Carlo simulation, we can 
account for all of the possible ns an electron with E e can 
produce and so predict the distribution of 77,5, which can 
then be compared with the experimentally determined 
distribution. 

One issue that we must accurately address in our sim- 
ulation is the sub-population of product electrons that 
escape the chamber before they are completely stopped 
(e.g., electrons created near the chamber's edge and mov- 
ing at large angles relative to the chamber's axis). Be- 
cause the number of sparks created by an escaping elec- 
tron is quite sensitive to its direction of travel, to model 
this type of event accurately requires our simulation to 
explore a large parameter space. To minimize the influ- 
ence of escaping electrons on the outcome of our simula- 
tion, we restrict our data set in the following way so that 
this sub-population of electrons is as small as possible: 
First, we note from Eq. Q that the most energetic (~ 
50 MeV) product electron will require about 7 cm of alu- 
minum to be stopped. Since each plate is approximately 
1 cm thick, about 7 plates will be required to stop such 
an electron that is directed straight downward. Also, we 
find empirically that a spark is created only if the elec- 
tron's path is inclined up to roughly 30° relative to the 
chamber's axis. Thus, by considering muon decays that 
occur only with the cylindrical "fiducial volume" defined 
to be in upper half of chamber with radius Rfid equal to 
half the plate radius i?, most of the spark-producing elec- 
trons will interact with enough aluminum to be stopped 
before escaping the chamber. 

For our simulation, [15 because charged particles lose 
negligible energy when passing through the chamber's 
spark gaps, we assume each muon and resulting product 
electron are stopped within the interior of a plate. At 
each muon-decay site, the product electron is assumed to 
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FIG. 4: Geometry of simulation, (a) Initial position r and 
angles and cj) for product electron; and (b) length / traversed 
within aluminum for stopped and escaping electrons. 



emanate downward, located a radial distance r (< Rfid) 
from the chamber's axis and directed at polar angle 
relative to the chamber's (downward) axis and azimuthal 
angle <\> as in Fig. [4j If the decay occurs a (perpendicu- 
lar) distance z Q above the plate's bottom surface and the 
electron's path length for being stopped in aluminum is 
/, then the number of sparks produced is given by 
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where the floor function return the highest integer less 
than or equal to its argument. Also, if I > l esc , the 
electron will escape the chamber before being stopped, 
where 



-r cos (j) + \J r 2 cos 2 <p + (R 2 — r 2 ) t p 



suit 



t<p -\~ tn 



(6) 



Hence, the algorithm for our simulation is as follows: 
Assume a value for m^c 2 . For the product electron ema- 
nating from a muon-decay site, randomly assign z Q , r, 
and (j) to be within the range — £p, — i?/2, — max 
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and — 27r, respectively, and also randomly assign the 
electron to have an initial kinetic energy E e within the 
range of — m /i c 2 /2. Then, using Eq. Q and (J6|, we 
calculate I and l esc . If I > l esci the electron's path length 
in the chamber is set to I = l esc . Finally, the number of 
sparks produced by this electron is found from Eq. ([5| 
and, based on Eq. (J2|, the number of times this event is 
expected to occur is assigned the value 
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N sim (n s ) = Dr(m^c z E e y(3 
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where D is a normalization constant. The factor of r 
is included because the differential area at radius r is 
dA = 2nrdr and so the number of muon decays at dif- 
fering radii is expected to be proportional to a geometric 
scaling factor of r. After iterating this procedure N*?^ 
times (typically = 10000), the accumulated N(ns) 

values are normalized so that J2 n N sirn (ns) = Nj;^, 
where N^ p is the total number of muon decays ob- 
served experimentally. Then, the predicted distribution 
N S im(ns) is compared with the experimentally deter- 
mined N exp (ns) as shown in Fig. [3] 

Finally, we iterate our simulation several times, where 
we assume a different value for m M c 2 each iteration. 
At the conclusion of each iteration, we first compare 
the experimentally determined N exp (ns) with the pre- 
dicted N sim (ns) by calculating the residuals r(ns) = 
N exp (ns) — N S i m (ns), and then determine the sum S 
of the squared residuals 



N e xp{ns) - N sim (n s ) 



(8) 



After concluding all of the iterations, we plot S vs. m M c 2 , 
and then obtain our "least-squares best-fit" value for the 
muon mass from the minimum of this plot. As shown 



in Fig. 
in Fig. 



the best-fit value derived from the data given 
is m M c 2 = 100 MeV. A weighted least-squares 



calculation yields the same result. 



VI. CONCLUSION 

As shown in Fig. [3J when our simulation is run under 
the assumption that m^c 2 = 100 MeV and max = 30°, 
N eX rp(ns) and N S i m (ns) are in very good agreement. To 
demonstrate the sensitivity of our approach to the as- 
sumed value of the muon mass, Fig. [6] displays the result 
of simply changing the value of ra^crused in the simula- 
tion to 90 MeV. We see that N exp (ns) and N 3 im(ris) now 
significantly disagree. Alternatively, if the plot of S vs. 
m M c 2 is generated under the assumption of max = 40°, 
the plot's minimum occurs at m^c 2 =105 MeV. Through 
such investigations as well as from noting that the min- 
imum region in Fig. [5] is relatively flat over the region 
from 95 to 105 MeV, we estimate our experimental value 
for m M c 2 to be (100±5) MeV, a range consistent with the 
accepted muon mass value of 106 MeV. 
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FIG. 5: Sum S of squared residuals vs. muon mass m M c 2 , 
assuming o-max 

— 30°. Minimum occurs are 100 MeV. 
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6: Comparison of experimental data (N*°p — 43) with 
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VII. APPENDIX: CONSTRUCTION DETAILS 

Spark Chamber: The design of the spark chamber is 
an adaptation of a spark chamber constructed by A.M. 
Sachs in the mid-1960s. [T6] The chamber is constructed 
of 3/8" thick, round aluminum plates. The chamber con- 
sists of 20 gaps and 21 plates. The 19 inner plates are 6" 
in diameter; the top and bottom plates are 7" in diame- 
ter. The conductive aluminum plates of the chamber are 
separated by 6" outer diameter, 5 1/4" inner diameter 
plexiglass spacer rings, 1/4" in thickness. It is essen- 
tial that the chamber be completely leak tight, and O- 
rings provide a simple cost effective means of sealing the 
chamber. Therefore, the gas seals between the aluminum 
plates and plexiglass spacers are made with O-rings. The 
plates of the chamber have grooves to accommodate 1 / 8" 
thick, 5 1/2" inner diameter O-rings. The 19 inner plates 
have 1/4" holes to allow gas to flow freely in the cham- 
ber. The chamber is held together by four threaded rods; 
these rods are electrically shielded from the pulsed plates 
with 1/2" O.D., 1/4" I.D. plexiglass sleeves. 

The chamber is filled with pure neon at slightly over 
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FIG. 7: Circuit diagram for capacitor bank. 



atmospheric pressure (approximately 1.1 atm) to ensure 
that any leakage in the chamber will be outward, not in- 
ward. This prevents any atmospheric gas impurities from 
leaking into the chamber. Electronegative gas impurities 
such as oxygen have detrimental effects on the chamber's 
operation. 

Scintillation Detectors: Each scintillation detector 
consists of a light-tight 20 cm x 20 cm x 2 cm plas- 
tic scintillator (Saint-Gobain BC-408) coupled to a blue- 
light sensitive photomultiplier tube (Hamamatsu R7400) . 
The detectors were purchased as assembled units (Saint- 
Gobain 8X8.8BC408/.5L-X). The PMTs of all three de- 
tectors are connected in parallel to a high-voltage sup- 
ply (Stanford Research Systems PS310) and operated at 
1000 V. The signal output of each PMT is passed through 
a fast (1 ns rise time) noninverting preamplifier (Ortec 
VT120C) with fixed 20 x gain. 

Coincidence/ Ant i- Coincidence Electronics: Each of 
the amplified PMT signals from the X, Y , and Z scin- 
tillation detectors is sent to an input of a multi-channel 
constant fraction discriminator (CFD) with shaping de- 
lay set to 4 ns. To avoid incorrect transient logic states, 
cable delays from detector to CFD are arranged to en- 
sure that Z arrives to the CFD before X and Y. We use 
the quad 200 MHz Ortec 935, a NIM-bin module whose 
negative-NIM logic output pulse widths are adjustable 
(typical value in our experiment is 300 ns). The X, Y", 
and Z NIM logic output signals from the CFD are input 
to each of two separate channels of a quad 4-Input Logic 
Unit (Ortec CO4020), where each channel is configured 
to perform the logical operation X • Y • Z. When this 



condition is met, each channel outputs a TTL pulse with 
width of 800 ns, one of which triggers the high-speed 
driver, while the other triggers the frame grabber board. 

High- Volt age Pulsing Electronics: Triggerable high- 
voltage pulses are applied to the spark chamber plates 
using the circuit shown in Fig. [7| Inside a homemade 
"capacitor bank," several thousand (typically 6000) Volts 
from a 10 kV high- voltage supply (Bertan 230- 10R) are 
applied to the anode of a hydrogen thyratron (Perkin 
Elmer HY-6) as well as the "high-voltage" side of each 
C = 500 pF capacitor (Sprague 20 DK-T5) in a parallel 
array of capacitors. The "low- voltage" side of each capac- 
itor is connected to a unique plate of the spark chamber 
as well as a current limiting R — 3 MQ resistor, whose op- 
posite end is grounded. This circuit is triggered by a com- 
mercially available Thyratron Driver (Perkin Elmer TM- 
27, modified to trigger on 1.5 V).[17 When the Driver 
receives a TTL pulse from the Logic Unit, it outputs a 2- 
/is wide pulse of amplitude 800 V with a rise time of less 
than 150 ns. This trigger pulse is passed into the home- 
made pulser box, where it is applied to the grid of the hy- 
drogen thyratron, driving the thyratron into conduction 
for 2 /is and causing the "low- voltage" sides of the capac- 
itors to drop to large negative voltage for a time period 
of approximately r = RC = (3 M^)(500 pF) = 1.5ms. 

Image Acquisition System: A monochrome progres- 
sive scan CCD camera head (Sony XC-55) with zoom 
lens (Navitar Zoom 7010) or else a camcorder (Canon 
ZR 40) with its autofocus option turned off is focused on 
the spark chamber directly as well as its perpendicular 
reflection in a mirror angled at 45°. With the camera's 
shutter speed at 1/30 s, the analog video signal is sent to 
a frame grabber board (National Instruments PCI-1405) 
plugged into a PCI expansion slot of a PC. The TTL 
output from the Logic Unit is connected the frame grab- 
bers TRIGGER input. Under the control of a Lab VIEW 
based program, when each TRIGGER pulse is received, 
the current frame acquired by the camera is stored under 
a unique filename as a jpeg file on the PC's hard drive. 



Acknowledgments 

We would like to thank Greg Eibel for invaluable as- 
sistance in machining our chamber and to Jim Brau and 
David Griffiths for helpful discussions. 



[1] Adrian C. Melissinos and Jim Napolitano, Experiments 
in Modern Physics, Second Edition (Academic Press, San 
Diego, 2003) pp. 399-409. 

[2] R. E. Hall, D. A. Lind, and R. A. Ristinen, "A simplified 
muon lifetime experiment for the instructional labora- 
tory," Am. J. Phys. 38, 1196-1200 (1970). 

[3] A. Owens and A. E. Macgregor, "Simple technique for 



determining the mean lifetime of the cosmic ray \i me- 
son," Am. J. Phys. 46, 859-860 (1978). 

[4] Roger J. Lewis, "Automatic measurement of the mean 
lifetime of the muon," Am. J. Phys. 50, 894-895 (1982). 

[5] T. Ward, M. Barker, J. Breeden, K. Komisarcik, M. 
Pickar, D. Wark, and J. Wiggins, "Laboratory study 
of the cosmic-ray lifetime," Am. J. Phys. 53, 542-546 



(1985). 

[6] Thomas Coan, Tiankuan Liu, and Jingbo Ye, "A com- 
pact apparatus for muon lifetime measurement and time 
dilation demonstration in the undergraduate laboratory," 
Am. J. Phys. 74, 161-164 (2006). 

[7] Nalini Easwar and Douglas A. Maclntire, "Study of the 
effect of relativistic time dilation on cosmic ray muon 
flux — An undergraduate modern physics experiment," 
Am. J. Phys. 59, 589-592 (1991). 

[8] C. Amsler, "The determination of the muon magnetic 
moment from cosmic rays," Am. J. Phys. 42, 1067-1069 
(1974). 

[9] "Development of Some New Experiments in Nuclear and 
Particle Physics," Am. J. Phys. 34, 9-18 (1966). 

[10] see Particle Detectors section of http://pdg.lbl.gov 

[11] David J. Griffiths, Introduction to Elementary Particles, 
Second Edition (Wiley- VCH, Berlin, 2008), p. 310-315. 

[12] W.R. Leo, Techniques for Nuclear and Particle Physics 
Experiments, Second Edition (Springer- Verlag, Berlin, 
1994) pp. 37-44. 

[13] L. Katz and A.S. Penfold, "Range-energy relations for 
electrons and the determination of beta-ray end-point en- 
ergies by absorption," Rev. Mod. Phys., 24, 28-44 (1 952) . 

[14] |http://pdg.lbl.gov/2009/AtomicNuclearProperties| 

[15] Here we present a simplified simulation algorithm, which 
can be duplicated by anyone with programming skill in 
some computer language. This algorithm is based on a 
"research- grade" simulation one of us [BB] wrote in Root, 
the standard analysis framework used in contemporary 
high-energy physics. We have found that the simplified 
simulation produces effectively the same results as our 
Root-based simulation. 

[16] A.M. Sachs, "Spark Chambers: A simplified system for 
the observation of particle trajectories in two types of 
chambers," Am. J. Phys. 35, 582-594 (1967). 

[17] Alternatively, North Star Research TT-G2 Thyratron 
Driver can be used. 



